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We study the phase behaviour of hard rods with length polydispersity, treated within a simplified 
version of the Onsager model. We give a detailed description of the unusual phase behaviour of the 
system when the rod length distribution has a "fat" (e.g. log-normal) tail up to some finite cutoff. 
The relatively large number of long rods in the system strongly influences the phase behaviour: 
the isotropic cloud curve, which defines the point where a nematic phase first occurs as density 
is increased, exhibits a kink; at this point the properties of the coexisting nematic shadow phase 
change discontinuously. A narrow three-phase isotropic-nematic-nematic coexistence region exists 
near the kink in the cloud curve, even though the length distribution is unimodal. A theoretical 
derivation of the isotropic cloud curve and nematic shadow curve, in the limit of large cutoff', is also 
given. The two curves are shown to collapse onto each other in the limit. The coexisting isotropic 
and nematic phases are essentially identical, the only difference being that the nematic contains a 
larger number of the longest rods; the longer rods are also the only ones that show any significant 
nematic ordering. Numerical results for finite but large cutoff support the theoretical predictions 
for the asymptotic scaling of all quantities with the cutoff length. 
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I. INTRODUCTION 

Rod-like particles in suspension can undergo a phase 
transition between an isotropic (I) phase without orienta- 
tional order and a nematic (N) phase where rods are pref- 
erentially oriented along the so-called nematic axis. This 
transition has been observed experimentally for chem- 
ical and biological systems^"^. Theoretically, the two 
main approaches for analysing systems of rod-like par- 
ticles are due to Maier and Saupe^ and Onsager''. The 
Maier-Saupe theory is based on a long range attractive 
interaction between particles (originally conceived of as 
due to van der Waals forces) and neglects density varia- 
tions. Such a theory is therefore appropriate for describ- 
ing thermotropic liquid crystals, in which orientational 
phase transitions are induced by changes in temperature. 
The Onsager theory, on the other hand, only considers 
hard core (short range) interactions between particles. 
Temperature then becomes just a trivial factor setting 
the energy scale, and the Onsager theory therefore de- 
scribes lyotropics: materials in which the phase transition 
is driven by a change in density of the system (at fixed 
temperature). The I-N phase transition is due, within 
the Onsager theory, to a competition between the orien- 
tational entropy — the tendency of rods to stay orienta- 
tionally disordered — and the packing entropy; the latter 
is due to the excluded volume interaction, and higher 
for aligned particles. The key simplification of Onsager 
theory is to consider the "Onsager limit" of thin rods, 
where the ratio D/L of the diameter and the length of 
the rods tends to zero. The free energy of the model can 
then be derived from a virial expansion truncated after 
the first nontrivial term, as all higher order terms turn 
out to be smaller by positive powers oi D/L. To find the 
actual free energy for a given rod number density, one 
needs to minimize over the orientational distribution of 
the rods; once this is done, the TN phase transition can 



be obtained by a standard double-tangent construction. 
Onsager solved the minimization problem by assuming a 
parametric form for the orientational distribution func- 
tion; even with only one variational parameter, he ob- 
tained results for the densities of the coexisting phases'"'* 
which are in good agreement with the numerically exact 
solution obtained later^. 

A comparison between experimental systems and the 
theoretical predictions of Onsager theory can be com- 
plicated due to 6.17. non-hard interactions or particles 
that are not perfectly rigid'^-'''^"-^^. A further impor- 
tant factor is polydispersity, i.e. a spread in particle 
lengths and/or diameters"^^'^, and we focus in this pa- 
per on the effects that length polydispersity has on the 
phase behaviour of hard rods. Polydispersity makes the 
analysis of the problem more difficult, but also leads 
to much richer phase behaviour^''. Even in simple bl- 
and tri-dispcrsc systems, i.e. mixtures of rods of two 
or three different lengths, phase separation into two ne- 
matic (N-N) phases and three-phase TN-N coexistence 
have been predicted theoretically^^'®'^^"^^, and observed 
experimentally^^. More generally, a pronounced broad- 
ening of the coexistence region and fractionation of the 
longer rods into the nematic phases are typical effects of 
length polydispersity^^'^^"^^'^. These effects are also seen 
in simplified models, such as the polydisperse Zwanzig 
modcl^*^'^^ — where rods are only allowed to point along 
one of three orthogonal directions — or the lattice model 
developed by Flory et al.^^~^^. Within the context of 
the lattice model, exponential and Poisson^^'^^ as well as 
Gaussian^"^ length distributions were analysed, with frac- 
tionation and a broadening of the coexistence region ob- 
served in all cases. The bidisperse case was also studied^"* 
and, despite the rather different theoretical approach, 
showed phase behaviour qualitatively similar to that pre- 
dicted within bidisperse Onsager theory^*'^^. 

It is clear from the above discussion that more work 
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is needed to understand the effects of length polydisper- 
sity on the phase behaviour of hard rods, especially with 
regard to the occurrence of the more "exotic" N-N and 
I-N-N phase coexistences. We recently began to inves- 
tigate these questions, using the "7^2 Onsager model" 
This is obtained from the conventional Onsager theory 
by simplifying the angular dependence of the excluded 
volume term in the free energy, truncating a series ex- 
pansion in Legendre polynomials after the first nontriv- 
ial term, which contains the second order Legendre poly- 
nomial V2 (cos 9) . The simplicity of the resulting model 
enabled us to investigate in detail the phase behaviour re- 
sulting from different choices of rod length distributions. 
For a Schulz distribution — which is unimodal, i.e. has a 
single peak — fractionation and broadening of the coex- 
istence region were found as expected, but neither N-N 
nor I-N-N coexistence appeared. These more complex 
features did occur, however, for bidisperse and bimodal 
distributions; the latter were modelled by a mixture of 
two Schulz distributions peaked at different lengths. For 
the bidisperse case the appearance of these features is 
in encouraging qualitative agreement with the results of 
Onsager theory, suggesting that the approximation made 
in constructing the V2 model may not be crucial. A de- 
tailed analysis of the phase diagram of the P2 Onsager 
model showed that N-N and I-N-N regions only appeared 
if the ratio of the rod lengths was sufhciently large while 
the fraction of long rods in the system remained small. 
(For the bimodal case the individual peaks in the dis- 
tribution also had to remain relatively narrow.) This 
suggested to us that N-N and I-N-N coexistence could 
also occur in unimodal length distributions that contain 
a larger number of long rods than the Schulz distribu- 
tion with its exponentially decaying tail. We therefore 
investigate, in the present paper, the behaviour of the V2 
Onsager model for rod length distributions with fat tails, 
i.e. decaying less than exponentially for large rod lengths. 
Our primary example will be the log-normal distribution; 
like the Schulz distribution, this has most of its weight 
around the average rod length, but contains a compar- 
atively larger number of much longer rods; it has also 
yielded interesting results in previous work on length- 
polydisperse homopolymers'^^''^^. Our main result is that 
I-N-N phase coexistence does indeed occur, although the 
topology of the phase diagram is rather different from 
the bidisperse case: the I-N-N region is very narrow in 
density and remains fully within the isotropic-nematic 
region, without being bordered by an N-N coexistence 
region. 

The present paper is structured as follows. In Sec. II 
we give a brief description of the polydisperse V2 On- 
sager model and of the phase coexistence equations. We 
also motivate why the presence of fat tails should cause 
unusual phase behaviour. To make the phase coexistence 
problem well-defined, a cutoff on the rod length distribu- 
tion needs to be introduced; this is of course also reason- 
able physically. Sec. Ill contains our numerical results for 
the phase diagram at finite cutoff, and we describe care- 



fully the unusual effects on phase behaviour caused by the 
presence of the long rods; it turns out that such effects are 
largely confined to a region of low densities where phase 
ordering would not yet be observed in systems with "con- 
ventional" length distributions. The second major part 
of the paper contains a theoretical treatment of the on- 
set of ncmatic ordering, in the limit of large length cutoff 
(Sec. IV). This reveals many unusual features; in the 
limit of infinite cutoff, for example, both the densities 
and the rod volume fractions of the isotropic cloud phase 
and the coexisting nematic shadow phase coincide. The 
theoretical predictions are compared to the results of our 
numerical calculations, and show good agreement. We 
conclude in Sec. V with a summary and outlook towards 
future work. Technical material is relegated to two ap- 
pendices. 

II. THE POLYDISPERSE T2 ONSAGER MODEL 

The 7^2 Onsager model is an approximate version of 
the Onsager model of hard rods. As in the Onsager the- 
ory, particles are modelled as rigid spherocylinders with 
hard core interaction. Phase transitions are driven by 
density rather than temperature; the latter simply fixes 
the energy scale and can be set to unity^'^^'^°'^°. 

We allow for length polydispersity, with the rod lengths 
L distributed according to a length distribution P{L), 
while assuming that all rod diameters D are equal. To 
be able to take the Onsager limit of thin rods simulta- 
neously for all rod lengths in P{L), wc introduce a ref- 
erence length Lq and consider the limit D/Lq — > while 
keeping the normalized lengths / = L/Lq constant. The 
thermodynamic state of the system is specified by the 
density distribution p{l, 0), where p{l, 9) dl (KI/(A-k) gives 
the number density of rods with (normalized) lengths 
in the range I . . .1 + dl and orientations in a solid angle 
dn around any direction at an angle with the ncmatic 
g^^jg30,20 then decompose p{l,9) into 

Pii,e)^pmie)^p„pm{9) (1) 

where we have isolated the overall density po and the 
normalized length distribution P{1). The angular prob- 
ability distributions Pi (9) are normalized in such a way 
that 

jdBPi{9)^l (2) 
using the shorthand 

d9 = —dcos9 
2 

and the convention that all angular integrations are over 
the range . . . tt. Conversely, the density pa is obtained 
by integrating p{l,9) over all I and 9, 

po= [dlde p{l, 9)^ [dl p{l) 
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The excess free energy of the model arises from the ex- 
cluded volume interaction of the rods. In the Onsager 
limit, the virial expansion can be truncated after the 
second-order term. The simplification of the Onsager 
model consists in simplifying the angular dependence of 
the resulting excluded volume term, by expanding in Leg- 
endre polynomials^ and truncating after the first non- 
trivial term. If all densities arc measured in units of 
l/[(7r/4)Z?i§], the excess free energy density then be- 
comes simply^^ 



7 ci 2 C2 2 



(3) 



where ci — 2 and C2 = 5/4. We have also defined the two 
moments of the density distribution 

pi = [dlde lp{l, 61) = I dl lp{l) I de Pi{6) (4) 



and 

P2 = JdldB lV2{cose) = j dl lp{l) j de 7^2 (cos 6*) Pi{9) 

(5) 

The first of these is just the rescaled rod volume fraction, 
pi = {L{)/D)(j}, while p2 contains information on the ori- 
entational order in the system. It is clear from Eq. (5) 
that for an isotropic phase {Pi{0) = 1) p2 vanishes, while 
< P2 < Pi for a nematic phase, whose angular distri- 
bution will be peaked around 9 = and 9 = tt; p2 '\s 
therefore a natural orientational order parameter for the 
system. In fact one can rewrite p2 as 



dl d9 lV2{cos9)p{l)Pi(9) = / dl lpil)S{l) (6) 



P2 



if one defines the order parameter S{1) as the average of 
the second Legendre polynomial for rods of given length. 



S{1)^ d9 V2{cos 9) Pi (9) 



(7) 



Eq. (6) can also be read as p2 = pq{IS{1)), where the 
average is taken over the normalized length distribution 
P(/) = p{l)/po. ^ 

As usual the ideal part of the free energy, comprising 
the ideal gas term and the entropy of mixing, is 



/id= d9dl pil,9)[lnp{l,9)-l] 



(8) 



Using the decomposition (1) and the normalization (2) 
we can then write the total free energy as 



A. Phase coexistence equations 

Since the rod orientations - as opposed to their lengths 
- are not conserved, the orientational probability distri- 
butions Pi{9) for a given density distribution p(l) over 
rod lengths have to be found by minimizing the free en- 
ergy (9), subject to the constraints (2). Introducing ap- 
propriate Lagrange multipliers we obtain the condi- 
tion 



SPm yf + jdl^il)jd9Pm)- 
p{l) [\nPi{9) + 1] + lp{l) [cipi - C2P2V2] + k{1) = 



Solving for Pi{9) gives 



exp(/c2P2'P2) 
/ d9 exp(/c2/027^2) 



(10) 



where we have introduced the shorthand 7^2 = P2 (cos 9) . 
This, together with (5) gives a self-consistency equation 



P2 



dl lp{l) 



Jd9 r2CXp{lc2P2V2] 
Jd9 CX.p{lc2P2'P2) 



(11^ 



which can be solved for p2. 

To calculate phase coexistences, we need the expres- 
sion for the chemical potentials and for the osmotic pres- 
sure. The chemical potentials are obtained by taking a 
derivative of the free energy with respect to the density 
distribution p(/). 



M0 = 



Sp{l) 

lnp(0 - In 



d9 cxp{-lcipi + IC2P2V2) (12) 



where we do not need to differentiate explicitly with re- 
spect to the Pi{9) since they are chosen to minimize 
the free energy. The osmotic pressure, derived from the 
Gibbs-Duhem relation, is 



n : 



J+ dl p{l)p{l) = po + 



ci 2 
^Pi 



C2 2 
JP2 



(13) 



Equality of the chemical potentials in a set of coexisting 
phases labelled by a = 1 ... P then leads to the following 
expression for the length distribution 



p^^nO = R{1) J d9 cxp(-;cipi"^ + IC2P'2'V2) (14) 

where R{1) is a function common to all phases. It can be 
j obtained by imposing the lever rule or particle number 



f^jdl p{l) [\np{l) - 1] + j dl p{l) j d9 Pi{9)\nPi{9) + 

(9) 

with / given by Eq. (3). 



conservation, 



Y^v^^)pi'^)ii) = p^^ii) 



(15) 
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where p^^\l) is the overall or 'parent' density distribution 
of the system and the fraction of the system volume 

occupied by phase a. This gives 



R{1) 



(16) 



Using Eqs. (10,14), the full density distributions over 
lengths and orientations are therefore 



p(o)(0exp(-/cipr+;c2pr^2) 



(a) 



J:,v('>) J de' exp(-/cipi' 



(fc) 



lC2P^^V'^, 



(17) 



We thus have a closed system of equations whose solu- 
tions determine the phase behaviour of the V2 Onsager 
model. For each of the P phases we have 3 unknowns, 
p'f'\ p^2^ and v^°'\ giving 3P unknowns in total. For each 
phase, p^"^ and p'^2^ obey the equations obtained by sub- 
stituting Eq. (17) into Eqs. (4,5); for /ci2"^ this leads back 
to Eq. (11). The remaining P equations are given by 
the equality of the osmotic pressure in all phases {P — 1 
equations) and the normalization condition u^"' ~ 1. 
The lever rule and the equality of chemical potentials are 
already ensured by having density distributions in the 
different phases of the form (17). 

So far everything holds for arbitrary parent distribu- 
tions p^^\l). We will denote by pg^"* the overall parent 

number density, so that P^°\l) = p'^^Hl) / Pq''^ is the nor- 
malized parent rod length distribution. We always as- 
sume that the average rod length in the parent is one; 
a different value could be absorbed into the reference 
length Lq. 



B. Fat-tailed distributions 

We now motivate why the phase behaviour of systems 
with rod length distributions with fat (less than expo- 
nentially decaying) tails should be unusual. To do so, it 
is useful to focus on the onset of the isotropic-nematic co- 
existence (the so called "isotropic cloud point"). Equal- 
ity of the chemical potentials (12) between an isotropic 
(p^(/)) and a nematic (p^(0) phase gives 

p^(0 = p\l)e''' [ M cxp [c2P^l{V2 - 1)] (18) 



where we have defined 



(19) 



and used p\ = 0. At the isotropic cloud point, i.e. at 
the onset of the phase coexistence between the isotropic 
parent and an infinitesimal amount of nematic "shadow" 
phase, p^{l) ~ Moreover, (3 should be positive 

since the nematic phase will have the larger rod volume 



fraction, pf > p\. The angular integral in Eq. (18), fi- 
nally, is bounded by unity since 7^2 — 1 < 0, and a Gaus- 
sian approximation around 9 = that is valid for large 
I shows that it varies only as a power law with I. Taken 
together, these facts imply that as soon as the parent dis- 
tribution p'"' (l) is less than exponentially decaying with 
I, the nematic length distribution would diverge for large 
I. In order to ensure convergent integrals for the density 
Pq and the rod volume fraction p^ we will therefore need 
to impose a maximum length cutoff Im on the parent dis- 
tribution. The need for such a cutoff suggests — and we 
will find this confirmed below — that even though there 
are only a very small number of long rods (since the in- 
tegral / dl p^°\l) converges) they can dominate the onset 
of phase coexistence from the isotropic side. Notice that, 
coming from the nematic, no such effects are expected. 
At the nematic cloud point, p^{l) = p'^^^l). Inverting 
Eq. (18) one sees then that the length distribution in the 
isotropic shadow phase is well behaved even for Im 00, 
as (3 is still positive and the angular integral is harmless. 

Motivated by these insights, we consider fat-tailed par- 
ent distributions in the rest of this paper. We expect, 
and indeed find, that the presence of the fat tail, and 
the value of the cutoff Im , will have a strong influence on 
the isotropic cloud point and the corresponding nematic 
shadow, while leaving the nematic cloud point essentially 
unaffected. Although many of our results are valid for 
general fat-tailed distributions, we will normally focus 
on log-normal distributions for definiteness. 



III. NUMERICAL RESULTS 

In this section we describe the unusual phase behaviour 
of systems with a log-normal parent length distribution 



p("'(0 



1 



\/2ttw^I 



exp 



(In^-^)^ 
2^ 



(20) 



with a finite cutoff, I <lm- The parameters p and w are 
chosen such that the average rod length is {I) = 1 and the 
polydispersity tr, i.e. the normalized standard deviation 
of distribution of rod lengths defined through 



(0^ 



(21) 



has the desired value; explicitly, this gives = ln(l+cr^) 
and p = —w^ 12. Note that for finite cutoff Z^, P^^^l) as 
given will not be normalized precisely to one, and (l) and 
(P) will differ slightly from the desired values of 1 and 
1 + a^. However, even for modest cutoffs the deviations 
arc very small. For instance, at cutoff Im = 50 and a = 
0.5, the integrals / dl P^°\l)r differ from their Z„ -> 00 
values 1, 1 and 1 + by values of order 10~^^, 10^^^, 
and 10~^^. In the following we will neglect these small 
corrections. 

Due to the numerical difficulties introduced by poly- 
dispersity, and compounded by the presence of a fat 
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tail in the parent distribution, we used the moment 
mcthod'^^"'^^ to solve numerically an approximate system 
of phase coexistence equations, applying the recently de- 
veloped adaptive method^"''^" to keep deviations from the 
exact solution as small as possible. We found that con- 
trolling such deviations was rather more difficult than for 
more "well-behaved" rod length distributions, and there- 
fore decided to integrate the moment method with the 
exact solution. In this mixed method, we effectively use 
the adaptive moment method to produce a good starting 
point for the numerical solution of the exact phase coexis- 
tence conditions. This method allows us to calculate the 
exact phase behaviour, while avoiding the convergence 
problems that result from a direct numerical attack on 
the exact phase coexistence conditions. 

A. Overview of the limiting theory 

Before showing the numerical results at finite cutoff 
Im, it is useful to anticipate some of the analytical re- 
sults which we will derive later on, in Sec. IV, for the 
large cutoff limit. This helps to understand and physi- 
cally interpret some of the unusual features of the phase 
diagram obtained at finite cutoff. 

The main feature affecting the physics of the system 
is obviously the presence of a large (relative to more 
strongly decaying distributions, e.g. of Schulz type) num- 
ber of long rods. Because of the weighting with / in the 
moment densities (4,5) appearing in the excess free en- 
ergy, it is plausible that these long rods can drive the 
phase ordering. Indeed, our theory will show that rods 
with / of the order of the cutoff become strongly or- 
dered already when the nematic phase first appears, and 
that this ordering drives the onset of phase coexistence 
to lower densities. The isotropic cloud point therefore 
depends on the cutoff, moving to lower densities as Im 
increases. However, the ordering of the longest rods is 
not sufficient to precipitate general order in the nematic 
phase. In fact, we find that at the relatively low isotropic 
cloud point densities caused by the presence of the long 
rods, the short rods with I of order unity are not yet 
able to order, staying completely disordered in the limit 
Im oo. The system thus has an unusual initial phase 
separation, in which the nematic phase is distinguish- 
able from the isotropic phase only because it contains a 
somewhat larger fraction of the longest rods, which are 
also nematically ordered; for the short rods, the den- 
sity distribution is essentially identical in both phases. 
At finite Im the contribution of the longest rods causes 
the nematic phase to have a larger rod volume fraction 
and average rod length than the isotropic phase. As Im 
increases, however, the differences turn out to vanish, 
and the number density and rod volume fraction of the 
nematic shadow phase become identical to those of the 
isotropic cloud phase. In a plot of these quantities ver- 
sus the polydispersity cr, the resulting cloud and shadow 



curves therefore coincide for Im oo. In fact, we will 
obtain this limiting curve explicitly as 

_ 1 1 

P'-P'- 4 (^^ 

Our theory also predicts the leading terms in the ap- 
proach to this asymptotic behaviour as Im increases. The 
number density and volume fraction of the isotropic and 
the nematic phase turn out to differ, in the limit, by 
terms proportional to ; here is the value of the mo- 
ment density (5) in the nematic phase, which for a log- 

1 /2 

normal parent scales as (lnZm)/Zm • The parameter (3 in 
Eq. (18) turns out to scale like (p^) to leading order. 
This decrease of (3 to zero allows the theory to remain 
well-defined, with finite values for the number density 
and rod volume fraction of the nematic for — > oo. In 
fact, the scaling of /3 is such that the two factors exp{f3l) 
and p^(l) = /9^°-*(0 in Eq. (18) approximately cancel for 
I = Im, with their product varying only as a power law 
in Im- In our numerical results at finite Im, this feature 
will be visible as a peak at I ~ Im in the density dis- 
tribution of the nematic shadow phase. The power law 
scaling of the weight of the peak is such that its con- 
tribution to the moment densities Pq and pf becomes 
negligible for Im — >■ oo, as pointed out above. On the 
other hand, higher moments of the nematic density dis- 
tribution, J dl p^{l)l"' with n>2, are dominated by the 
contribution from the peak and actually diverge with Im- 
Similar unusual features also appear in, for example, the 
Flory-Huggins theory of polymers with log-normal chain 
length distribution'^^''^^. 

B. Overall phase diagram topology 

In Fig. 1 we show the phase diagram obtained for a log- 
normal parent distribution with cutoff Im ~ 100. Plotted 
is the density of the parent at which phase transitions 
occur; the lines separating one- and two-phase regions 
are therefore cloud curves. 

Vertically, the phase diagram is broadly divided into 
two regions: in the first one, at low polydispersity, the 
phase behaviour mirrors closely that of a system with 
a Schulz distribution of rod lengths^°, exhibiting an I- 
N coexistence region in the middle and single-phase I 
and N regions to the left and right. The isotropic cloud 
curve essentially coincides, at these small values of a, 
with the one found for a Schulz length distribution^'' and 
indicated by the dashed line in Fig. 1 . At the other end of 
the I-N coexistence region the nematic cloud curve (not 
shown) is, as expected from the discussion in Sec. II B, 
essentially unaffected by the presence of long rods and 
remains close to its analogue for a Schulz distribution in 
the whole range of a. 

The second region of the phase diagram, at higher a, 
is more strongly affected by the relatively large number 
of long rods in the log-normal parent distribution, and 
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FIG. 1. Phase diagram for a log- normal parent with cut- 
off = 100. Shown is the polydispersity cr plotted against 
the density of the parent Pq"' at which phase transitions oc- 
cur. From left to right we find a single-phase isotropic phase, 
the I-N coexistence region and then the single-phase nematic; 
the coexistence region is bounded from left and right by the 
isotropic and nematic cloud curves, respectively. Inside the 
I-N coexistence region we observe a narrow three-phase I-N-N 
region. This is shown in more detail in the inset and lies close 
to the isotropic cloud curve (dashed line) that would result 
for a Schulz distribution. Above the three-phase region one 
has, instead of N-N phase separation, a crossover from a cut- 
off-dependent regime (located to the right of the isotropic 
cloud curve) to a regime in which the phase behaviour is al- 
most unaffected by Im ; the location of the crossover is close to 
the dashed line, i.e. the "conventional" (Schulz) cloud curve. 



therefore also by the value of the cutoff At the bound- 
ary between the two regions, at ct « 0.24, the cloud curve 
(see Fig.l) has a kink. Correspondingly, the properties of 
the coexisting nematic shadow phases change discontinu- 
ously at this point, from a "conventional" nematic phase 
similar to that found for a Schulz distribution to one 
containing more, and more strongly ordered, long rods; 
see Figs. 2 to 4 below. Intuitively, below the threshold 
value of a the parent distribution is narrow enough to 
keep the number of long rods small, to the point where 
phase behaviour is not significantly affected, while above 
the threshold the influence of these long rods can become 
dominant. 

Precisely at the kink in the isotropic cloud curve, the 
isotropic phase coexists with two different shadow ne- 
matic phases. As expected, we therefore find a three- 
phase region beginning at this point. Crossing the 
isotropic cloud curve from the left, the system thus first 
separates into an isotropic phase and an "unusual" ne- 
matic (Ni) with many long rods. Upon increasing den- 
sity, there is a narrow interval of three-phase I-N1-N2 co- 
existence, with N2 being a "conventional" nematic phase. 
Beyond this point, the system behaves in the usual way, 
showing an I-N2 coexistence and finally a single nematic 
phase. On increasing cr, however, the three-phase region 



eventually closes off and the strict transition between 
Ni and N2 is replaced by a pronounced crossover in the 
shape of the density distribution of the nematic phase; 
see Sec. Ill D below. The crossover line in the phase di- 
agram essentially continues the three-phase region, and 
both are close to the isotropic cloud curve for a Schulz 
distribution (dashed line in Fig. 1), which is representa- 
tive also of other parent distributions without fat tails. 
We can thus say that it is essentially only to the left of 
this line in the phase diagram that the fat tail of the log- 
normal parent, and therefore also the value of the cutoff, 
significantly affect phase behaviour; one of the resulting 
effects is of course that, above the threshold value of tr, 
the actual isotropic cloud curve is moved significantly to 
the left of its counterpart for a Schulz distribution. The 
unaffected region of the phase diagram includes, in par- 
ticular, the nematic cloud curve as expected from the 
discussion in Sec. II B. 

The fact that the three-phase region terminates at a fi- 
nite value of cr, where a strict phase coexistence involving 
two nematic phases turns into a crossover in the proper- 
tics of a single nematic phase, indicates that there must 
be a nematic-nematic critical point on the I-N-N phase 
boundary. However, because the three-phase region itself 
is extremely narrow, we have not been able to locate this 
point numerically. 

C. Cloud and shadow curves 

We now discuss in more detail the effect of the long 
rods on the isotropic cloud and nematic shadow curves; 
as shown above, the nematic cloud and isotropic shadow 
curves do not need to be considered further since they are 
essentially unaffected by the presence of the long rods. 
Fig. 2 shows the isotropic cloud (main plot) and nematic 
shadow (inset) curves, in the p^-a representation of plot- 
ting polydispersity versus the number densities po in the 
cloud and shadow phases. The kink in the cloud curves 
can clearly be seen, as can the associated discontinuity 
in the shadow curves. Upon increasing the cutoff length, 
the curves also move towards the predicted limit curve 
for /,„ — > 00. (We defer a discussion of the numerical re- 
sults for the actual rate of convergence to Sec. IV F.) The 
same is true if, as shown in Fig. 3, we use the rescaled 
volume fractions p\ rather than the number densities to 
represent the properties of the cloud and shadow phases. 
Interestingly, at the discontinuity in the shadow curves 
the volume fraction of the nematic phase becomes 
larger (Fig. 3) as cr is increased, while the number den- 
sity becomes much smaller and drops below that of its 
isotropic counterpart (Fig. 2). This implies an increase in 
the average rod length (/) = / of the nematic phase, 
i.e. a strong fractionation effect which is rather more pro- 
nounced than for a system with a Schulz distribution. 
This is shown explicitly in Fig. 4, where we plot the aver- 
age rod length in the nematic shadow phases against the 
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FIG. 2. Isotropic cloud curve and nematic shadow curve 
(in the inset) in the po-a representation, for the different val- 
ues of the cutoff Im shown in the legend. Even for the modest 
values of Im used, the evolution of the curves with increas- 
ing Im is compatible with the predicted convergence to the 
same limiting curve (thin solid line) for both isotropic cloud 
and nematic shadow. Note the discontinuity in each shadow 
curve, which corresponds to the kink in the associated cloud 
curve; here the nematic phase jumps from a "normal" phase 
to one dominated by the long rods. Below the kink or dis- 
continuity, the cloud and shadow curves are practically cut- 
off-independent and essentially identical to those found for a 
Schulz distribution. 



polydispersity. The discontinuities in the shadow curves 
for pi and po give a similar discontinuity here; below 
it, the results are very close to those for a Schulz dis- 
tribution, while for higher values of a the rod length in 
the nematic shadow phase is considerably larger, indi- 
cating that the nematic phase has become enriched in 
the long rods. (Note that the coexisting isotropic cloud 
phase always has average rod length one, since its density 
distribution is identical to that of the parent.) As a is in- 
creased, the fractionation first increases, but then reaches 
a maximum and decreases again. This can be understood 
by looking back at Figs. 2 and 3: the nematic shadow's 
number density p^ initially decreases more quickly with 
increasing a than its rod volume fraction p^ , leading to 
the observed increase in (1) = p^ /p^; but eventually Pq 
becomes almost independent of a and (l) then decreases 
with pf . It is intriguing that the maximum fractionation 
is reached at roughly the same value of cr at which the 
three-phase region disappears (compare Fig. 1), but we 
have no explanation for this at present. 

For completeness. Fig. 4 also shows the average rod 
length in the isotropic shadow; recall that this is defined 
as the isotropic phase which coexists with the nematic 
cloud phase at the upper limit (in density) of the phase 
separation region. As expected, the results show no pro- 
nounced effect from the presence of long rods and are 
close to those for a Schulz distribution. 

More insight into the properties of the isotropic cloud 
and nematic shadow phases can be gained by plotting 




PI = (Lo/D) 01 

FIG. 3. Analogue of Fig. 2, but in the pi-a representation, 
i. e. with (rescaled) volume fractions pi used instead of number 
densities po to indicate the properties of the phases. Note 
the large increase in p^ at the discontinuity of the nematic 
shadow curves, which contrasts with a decrease in pQ (see 
Fig. 2). Below the discontinuity, the results are again very 
similar to those for a Schulz distribution. 



their density distributions p{l). Fig. 5 shows exemplary 
results for = 100 and a value of a above the disconti- 
nuity in the shadow curve, i.e. in the region of the phase 
diagram where the long rods have a significant effect. As 
anticipated in Sec. Ill A, the density distribution of the 
nematic shadow for rod lengths of order unity has a shape 
similar to that of the isotropic cloud phase, while exhibit- 
ing an exponential increase for rod lengths / comparable 
to the cutoff Im- Note that the total density in this long 
rods part of the distribution is much smaller (by almost 
two orders of magnitude in the example) than the overall 
nematic density; but the fact that the lengths themselves 
are large makes this sufficient to give the large average 
rod lengths (l) that we found above. With increasing cut- 
off the peak moves to larger lengths while, as we will see 
later in more detail, its weight decreases. This is why, 
in the limit Im oo, the isotropic cloud and nematic 
shadow phase become indistinguishable in their number 
density and rod volume fraction, i.e. in their moments po 
and pi; as a consequence, the average rod length in the 
nematic shadow phase also tends to unity (the value in 
the parent) in the limit. 

Looking ahead to our theoretical treatment below, the 
fact that the nematic density distribution has two max- 
ima, one for / of order unity and a small second peak for 
I = Im, will allow us to split integrals over rod lengths 
into corresponding short and long rod parts. The log- 
scale representation of the density distribution in the ne- 
matic cloud phase (Fig. 5, bottom) supports the viability 
of this approach, showing a clear dip of p^{l) to negligible 
values between the short and long rod regimes. The fact 
that the decay from the second peak in p^{l): for I < Im, 
is close to exponential will further simplify matters, al- 
lowing us to replace non-exponential factors in the long 
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FIG. 4. The average rod length (l) in the nematic and 
isotropic shadow phases (solid) at cutoff Im = 100, compared 
with the same results for a system with a Schulz distribution 
(dashed) . The isotropic shadow, i. e. the phase which initially 
coexists with the nematic when the density is reduced from 
large values, is in both cases the one at lower (Z); as expected, 
it is essentially unaffected by the presence of long rods so that 
the results for the log-normal and the Schulz distributions are 
very similar. The same is true for the nematic shadow phase 
below the threshold value of a; above, the average rod length 
is much larger. The nematic phase here has an enhanced 
concentration of the longest rods, leading to pronounced frac- 
tionation. 



rod integrals by their values I = Im- 

Finally, we comment on the strength of orientational 
order in the nematic cloud phase. In Sec. Ill A we an- 
ticipated that in the limit of large cutoff this actually 
vanishes, as indicated by the convergence of to zero. 
This implies that rods with lengths I of order unity are, 
in the limit, orientationally disordered: Eq. (10) gives a 
uniform orientational distribution when p^l — > 0. More 
specifically, Eq. (10) implies that only rods with / of or- 
der l/p2 and greater show significant orientational or- 
der, and P2 for /,„ — > 00 implies that this "order- 
ing length" I/P2 diverges (while remaining <^ Im] see 
Sec. IV D). The convergence of p^ to zero is rather slow 
however, and for the modest values of /„i used so far 
is still of order unity. For example, the nematic shadow 
phase at CT = 0.5 and /,„ = 100 has « 2 and therefore 
noticeable orientational order even for rod lengths I of 
order one. 



D. Long rod effects in the coexistence region 

In Sec. IIIB, we reported that the region of the phase 
diagram most strongly affected by the long rods is lo- 
cated above the treshold value of a and between the ac- 
tual isotropic cloud curve and that which is obtained for 
distributions without fat tails, e.g. of Schulz type. We 
now explore this region in more detail, focussing on the 
regime of large a where there is no three-phase I-N-N co- 
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FIG. 5. Density distribution of the isotropic cloud phase 
(solid) and the nematic shadow phase (dashed), at the 
isotropic cloud point for Im = 100 and o — 0.63, in linear 
(top left and right) and log-log (bottom) scale. The top right 
panel shows the region of large I, where the nematic shadow 
distribution exhibits an exponentially increasing regime. The 
bottom panel demonstrates that, given the very small density 
of such long rods in the isotropic phase, their enrichment in 
the nematic phase is quite dramatic. 



existence but instead a pronounced change in the prop- 
erties of the nematic phase as the coexistence region is 
traversed. 

While at the isotropic cloud point it is only the longer 
rods in the nematic phase that exhibit noticeable orien- 
tational order (at least for large cutoff see end of 
previous section), we expect also the shorter rods to be- 
come ordered eventually as density is increased. As the 
effect of the longest rods on the nematic phase dimin- 
ishes, fractionation effects are then expected to become 
less pronounced, leading to a reduced average rod length 
in the nematic. A plot of the average rod lengths in the 
isotropic and nematic phases against the parent density 
across the coexistence region (Fig. 6) shows this effect 
clearly. At parent densities Pq""* « 2 we see a crossover, 
where the average rod length in the nematic drops to the 
values typical for phase separation from a parent distri- 
bution without a fat tail; as expected, the crossover den- 
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FIG. 6. Average rod length in the isotropic and ne- 
matic phase across the coexistence region at a — 0.7 and 
Im = 100. Notice the strong fractionation in the first, "long 
rod-dominated" , part of the coexistence region up to parent 
densities around p^^^ — 2, while above the behaviour resem- 
bles that for a parent with e.g. a Schulz distribution. Notice 
also the non-monotonicity of the average rod length in the 
nematic phase, which is consistent with our theory for the 
large-cutoff limit (see text). 




FIG. 7. Number density of the isotropic and nematic 
phases across the coexistence region for a = 0.7 and Im ~ 100, 
plotted against the parent density. Notice that in the long 
rod-dominated region below Pq"' « 2, the density pj, of 
the isotropic phase is almost identical to that of the parent 
(dashed line). In the same regime, the nematic density pQ 
shows a pronounced increase, before crossing over to the more 
gradual variation typical of systems without fat-tailed length 
distributions. 



sity is close to the isotropic cloud point (pg ~ 2.1) for a 
Schulz distribution with the same polydispersity. Notice 
that the average rod length in the nematic phase, {l)^ , is 
non-monotonic in parent density, initially showing an in- 
crease when moving away from the isotropic cloud point 
but then decreasing again. This is consistent with our 
theory for the large cutofF-limit, which predicts that for 
Im oo, (l)^ = {ly (= 1 with our convention for the 
parent's average rod length) at the isotropic cloud point. 
On increasing the density, (l)^ must then first increase as 
the system begins to fractionate, before reducing again 
to unity as the nematic cloud point is approached. 

The variation of the density of the two phases across 
the coexistence region (Fig. 7) similarly shows the 
crossover between two different regimes for the nematic 
phase: in the low-density regime, the density of the ne- 
matic increases rather quickly with the parent density, 
before crossing over to the more gradual variation that is 
familiar from systems with length distributions without 
fat tails. Interestingly, the density of the isotropic phase 
is essentially identical to that of the parent throughout 
the low-density region where the phase separation be- 
haviour is dominated by the long rods. Since the overall 
density of the system must equal that of the parent, the 
fractional system volume occupied by the nematic phase 
must then be very small. This is indeed what we find, 
and the representation in Fig. 8 of the phase behaviour 
provides probably the clearest visual demonstration of 
the crossover between the low-density regime dominated 
by long rods and the region of more conventional phase 
behaviour at higher densities. 



The most detailed information about the phase sepa- 
ration is obtained by analysing the density distributions 
in the coexisting phases. Focussing on the low-density 
part of the coexistence region that is dominated by the 
effect of the long rods, we find as expected from Figs. 7 
and 8 that the density distribution of the isotropic phase 
(not shown) is essentially identical to that of the parent, 
except for a reduction in the already small density of 
long rods. For the nematic phase, we saw in Fig. 5 that 
p^{l) has a peak a.t I = Im at the isotropic cloud point. 
When we move into the coexistence region by increas- 
ing the parent density (Fig. 9), the density distribution 
acquires a bimodal shape as the second peak gradually 
moves to smaller values of I. Eventually, at the point 
where the system crosses over to conventional phase be- 
haviour, this second peak merges with the main peak at 
I of order unity. Before this happens, the nematic phase 
remains significantly enriched, compared to the parent 
in some of the longer rods. Since the overall density 
distribution must be preserved in phase separation (see 
Eq. (15)), this again makes plausible that the fractional 
system volume occupied by nematic phase has to be very 
small in this regime: there simply are not enough long 
rods to allow formation of a significant amount of a ne- 
matic phase that contains a much higher proportion of 
such long rods. 

Finally, we may ask how the properties of the coexist- 
ing phases depend on the cutoff Im, at some fixed point 
in the region of the phase diagram where the long rods 
significantly affect the behaviour (i.e. to the left of the 
isotropic cloud curve for rod length distributions with- 
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FIG. 8. Fractions of system volume occupied by the 
isotropic (dashed) and nematic (sohd) phases, plotted against 
parent density across the coexistence region for a — 0.7 and 
l,n = 100. On a logarithmic scale (inset) it is clear the in the 
regime where the phase separation is dominated by the long 
rods, the fractional volume of the nematic is extremely small. 



out fat tails). Naively, one might expect that the cutoff 
vifould play an important role not just at the isotropic 
cloud point, but throughout this long rod-dominated part 
of the phase diagram. However, this is not so. It turns 
out, for example, that the peak in p^{l) at large / (but 
below the cutoff) docs not change significantly as l„i in- 
creases and eventually stabilizes. The same applies to all 
other properties of the coexisting phases, including the 
fractional volume of the nematic phase; the latter re- 
mains nonzero (though very small) as Im becomes large. 
Even while the isotropic cloud curve moves towards the 
limiting curve for infinite cutoff (Eq. (50)), widening the 
phase coexistence region, the phase separation behaviour 
in the region beyond the cloud point therefore eventually 
becomes independent of /,„. 

IV. THEORY FOR ISOTROPIC CLOUD POINT 

Having discussed the phase behaviour of systems with 
a log-normal parent distribution with a finite length cut- 
off Im, as obtained from our numerical calculations, we 
now turn to the limit of large Im where the onset of ne- 
matic ordering, i.e. the properties of the isotropic cloud 
and nematic shadow phases, are amenable to theoret- 
ical treatment. Even though any physical system will 
necessarily have a finite largest rod length and therefore 
a finite Im-, the limiting case Im oo is useful in un- 
derstanding some of the physics that we have described 
in the previous sections. In particular, we can extract 
a large amount of information about the characteristics 
of the coexisting phases, such as the scaling of all the 
densities of the coexisting phases with Z,„, which can be 
compared directly with our numerical results. The cloud 




FIG. 9. Plot of the density distribution of the nematic 
phase for short (left) and long (right) rods, for a range of 
parent densities across the coexistence region at cr = 0.28 and 
Im ~ 100. With increasing parent density, i.e. when moving 
away from the isotropic cloud point, the first (left) peak of 
the distribution grows and approaches the shape of the par- 
ent. The second peak at Im, on the other hand, first decreases 
in height and then gives way to a maximum at I < Im which 
moves towards smaller /. The crossover to the region in the 
phase diagram where long rods no longer significantly affect 
the phase behaviour takes place when this maximum merges 
with the main maximum on the left. 



and shadow curves can be deduced explicitly in the limit 
Im — oo and actually coincide. The moment density 
P2 of the nematic shadow phase vanishes in the limit. 
Since specifies the degree of orientational ordering, 
one may doubt whether this solution can justifiably be 
referred to as a nematic phase. However, recall from 
Eq. (6) that P2 is just Pq{S{1)1), where the average is 
taken over the normalized length distribution in the ne- 
matic, P^{1) = p^{1)/pq. a small value of then im- 
plies that S(l), the orientational order parameter for rods 
of length I defined in Eq. (7), must be close to zero wher- 
ever P^{1) is non- vanishing. In regions where P^{1) has 
little weight, on the other hand, the value of S{1) can be 
close to one. This is exactly what happens in the nematic 
shadow phase in the limit Im 00. For I of order unity 
P^{1) has nonzero weight as can be seen from the plot of 
= PoP^il) in Fig. 5. Correspondingly, S{1) ^ in 
this regime of rod lengths, implying that the short rods 
are orientationally disordered. For the longest rods with 
I of order Im, on the other hand, P^{1) has negligible 
weight for Im — *■ 00, and we will find that these rods are 
strongly ordered, with S{1) — > 1. The contributions to 
P2 from both regimes vanish in the limit, although non- 
trivially they are of the same order; the fact that the 
length distribution has much more weight on short rods 
is compensated by the smallness of S{1) there. 

In the next sections we will establish, by a process of 
elimination, the overall limiting behaviour of the relevant 
quantities as Im becomes large. Once we have derived 
the system of equations to be solved for finite Im, we 
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will decompose all integrals over / into two parts, one 
for lengths of order unity, I — . . .1 and one for the 
large lengths / = I . . Am- Here I is chosen such that 
?m S> Z 3> 1; formally, we first take l,n oo and then 
I ^ oo. The decomposition into the two parts allows us 
to simplify all integrals, making use of the fact that the 
short rods are orientationally disordered while the long 
ones are strongly ordered. From this we eventually derive 
the solution for the isotropic cloud point and the scaling 
of all quantities with Im- 

A. Asymptotic behaviour of (] and p2 

Let us start our analysis with a discussion of the be- 
haviour, in the limit 1^ oo, of the key quantity /3; from 
Eq. (18) it follows that its value has a strong effect on 
the properties of the nematic shadow phase. Explicitly, 
we have from Eq. (18) that the moments of the nematic 
shadow phase are 

^PojdlM P^°\l)e^'e'f^'>^ (22) 

p^ = PoJdld9 P^^'^Ole^' e'f^'^ (23) 

p^ = poJdldB P(°)(0 I 7^2 (cos 0) e'^C') (24) 

where we have defined f(6) as 

f{e) = C2P^{V2 - 1) (25) 

and have used the fact that the isotropic phase coincides 
with the parent distribution p^°^(0 at the isotropic cloud 
point; Po = Jdl p^^^l) is the parent density (we drop 

the superscript on Pq"' for brevity). An important con- 
straint on the dependence of /3 on 1^ is that all the above 
moment densities must remain finite as 1^ — > oo. This 
follows from the equality of the osmotic pressure (13), 
which gives 

Po + lcipl^p^ + \ci{p^f-^-C2{p^)' (26) 

On the l.h.s. we have used again that the isotropic cloud 
phase coincides with the parent, giving p^ — po and also 
Pi = Pq = Po since the average rod length in the parent 
is unity. All terms on the l.h.s. of Eq. (13) are finite. 
The r.h.s. is> p^ + (ci - C2)ip^y/2 since p^ < from 
Eqs. (23,24). Using ci > C2 it follows that both terms 
are positive, so that Pq and pf , and therefore also p^, 
must remain finite as claimed. As in Sec. II B, we can 
now deduce that a length cutoff Im must be imposed on 
all integrals over I whenever the parent distribution is 
less than exponentially decaying at large /. This follows 
because /3 should be positive in order to have a greater 
rod volume fraction in the nematic phase than in the 
isotropic, i.e. p^ > p\. Without a cutoff on I, such a 



positive value of f3 in Eqs. (22,23) would lead to diver- 
gent integrals for p^ and p^ , violating the constraint that 
these moment densities must be finite. (The angular in- 
tegrals cannot counteract the effect of the exponential 
factor exp{pi), since they vary at most as power laws 
with I.) The cutoff Im will be assumed large and later 
taken to infinity to study the limiting behaviour for large 
cutoffs. Initially, we develop the theory for general parent 
distributions with fat tails, but then eventually special- 
ize to the log-normal case in order to compare with our 
numerical results for finite cutoffs'^^. 

We are now in a position to narrow down the limiting 
behaviour of /3 as ^ oo. With f3 being positive, there 
are three possibilities: (3 — > oo, (3 — > const, and /3 ^ 0. 
The first case can be excluded from Eq. (19), which using 
p\ = Po becomes 

P = -ci {p^ - po) + C2p^ (27) 

We have already established that all terms on the r.h.s. of 
this equation remain finite, so that (3 cannot diverge. In 
the second case, /? — ^ const. > 0, the integral in Eq. (22) 
for p^ will diverge as l,n — > oo, giving p^/po oo and 
therefore po to keep p^ finite. But then the l.h.s. of 
the osmotic pressure equality (26) becomes ~ po, while 
the r.h.s. is > p^ -I- (ci — C2)(pf )^/2 > Pq as explained 
above; the pressure equality can therefore not be satis- 
fied since Pq ^ po. We therefore conclude that the only 
possible limit behaviour is P —f for Im ^ oo. This 
means that in the "short rods" region, i.e. for fixed val- 
ues of the rod length I that do not grow with Im, (31 is a 
small parameter in which we can expand. The limiting 
behaviour of the product [3lm is more difficult to deter- 
mine, although it is plausible that if this quantity is not 
large enough then the supposed nematic shadow phase 
will behave as if /3 = exactly and thus become fully 
identical to the isotropic cloud phase in the limit, imply- 
ing that no genuine phase coexistence has been found. 
We therefore assume in the following that j3lm ~^ oo for 
Im oo, and will find support for this in the numerical 
results presented in Sec. IV F. 

Let us now turn to the behaviour of p^ for Im oo. 
Having shown that it cannot diverge, the two possibili- 
ties are that p^ const. > or p^ ^ 0. In the first case 
we show in App. A that the problem can be drastically 
reduced to a single equation in p^. The only solution 
of this equation is p^ ~ 0, in contradicition with the 
assumed finite limit value. We therefore conclude that 
both (3 and p^ vanish in the limit of infinite cutoff. As 
explained above, we will see that the nematic shadow 
phase remains a bona fide nematic in spite of the fact 
that p^ — > 0: the longest rods exhibit strong nematic 
order but their contribution to p^ vanishes in the limit. 
As pointed out earlier, the fact that p^ also demon- 
strates that the short rods in the nematic shadow phase 
are orientationally disordered. 
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B. Setup of the equations 

We now proceed to simplify the equation for the 
isotropic cloud and nematic shadow phases, using that 
(i and both converge to zero for large Im while pim 
diverges. The equations to be solved are (22,23,24) for 
the moments of the nematic shadow phase, the osmotic 
pressure equality (26) and the definition (27) of /3. These 
five equations determine the five unknowns po, Po ' ^ 
P2 and /? as a function of the cutoff Im- Equivalently, we 
can use Eqs. (22,23) to eliminate pg and pf , solving the 
three remaining equations for po, and /3. 

We now proceed to implement our scheme of breaking 
the integrals over / into two intervals, namely I = . . .1 
and I = I . . .Im; where 1 ^ Z <C ^m. Calling /q the 
contribution of the long rods to p^ , Eq. (22) becomes in 
this way 



pN^p^y^ dl p(°)(0e'''y55eW)+Jo 

where 

Jo = po ^'"^ dl {l)e^' J dd e'-^(«) (28) 

In the integral over the short rods we can expand in the 
small parameters pi and p^l- Since from Eq. (25) f{0)l 
is proportional to p^l, an expansion to linear order gives 



dl P(°)(0 d0 {l + f3l + fid)l} + Io 



N 

Po Po 



Having isolated the contribution from the long rods, we 
can now extend the integration to oo in the convergent 
Z-integral in the first term; the error from this approxi- 
mation turns out to be negligible compared to the cutoff- 
dependent terms that we focus on. This gives, using that 
the average parental rod length is unity. 



Po = PO + W- C2P2)P0 + lo 



(29) 



The same procedure can be applied to Eqs. (23,24) and 
yields 



Pi =Po + (/3-C2P^')(Z')po + /: 

P2 = y^'0(' )P2 +h 

where the long rod parts are 

/i=po t dl /p("'(0e'3V-^W 



h=Poj^^ d//p('')(0e'3'7'2e'/(^) 



(30) 
(31) 



(32) 
(33) 



and the averages {I"') are taken over the normalized par- 
ent distribution 



{r) ^ dl rp^°\i) 



The long rod contributions /q, /i, I2 remain undeter- 
mined so far. They must certainly be finite for Im 00, 
because the moments are, but in fact will turn out to 
converge to zero. For I2 this is easy to see from Eq. (31), 
which implies that I2 < P2 ■ We can further exclude 
the possibility that I2 ^ p^, which would imply from 
Eq. (31) that 

{l^)po = 4 

It is easy to show that this is just the spinodal condi- 
tion for the isotropic phase in the P2 Onsager model, 
while the onset of isotropic-nematic phase coexistence 
must happen at a density below the spinodal instability. 
Thus, I2 must be of the same order as p^, which means 
that short and long rods make comparable contributions 
to P2 as anticipated above. 

To evaluate the integrals for Iq, Ii and I2 more quanti- 
tatively, we note that they all contain the exponentially 
diverging factor exp(/3Z) which dominates the behaviour 
of the integrands since (3l„i ^00. To leading order 
we can thus approximate the integrands P^^^ {l)w{l)e^'' , 
where w{l) is one of the weight functions in the above 
Eqs. (28,32,33), by p("' (/„)w(/,„)e'^' and integrate. For 
the log-normal parent distribution we show in App. B 
that this approximation indeed becomes exact for — > 
00. For more quickly decaying distributions, for example 
P^°^l) - exp(-/") with < a < 1, a slightly better 
approximation is needed (see App. B), but the integrals 
remain dominated by their values at I = Im- 

It follows from the above arguments that /q is smaller 
than Ii by a factor of l//m for large Im, and can therefore 
be neglected to leading order. Likewise, the ratio of I2 
and /i becomes simply 



I2 
h 



Jd9 P2e 



jde 



(34) 



which is in fact identical to S{lm), the orientational order 
parameter for the longest rods, and therefore lies between 
and 1. Using the definition of s, Eq. (31) can be written 
as 



p^ = ^Po(Z>^ + ./i 
5 



(35) 



which will be useful shortly. 

So far we have only dealt with the equations for the 
moment densities pg , pf and p^ of the nematic phase. 
We now add Eq. (27) for (3; inserting Eq. (30) and bring- 
ing all terms in (3 to the left gives 

P{l + ci {P)po) = C2p^ (1 + ci {P)po) - cih (36) 

This implies (3 < C2P2 , and since /3lm ^ 00 for Im ^ 00 
we also have p^lm — > 00. Recalhng f{6) = C2P2 {'P2 — 1), 
we therefore see from (34) that the long rods have strong 
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orientational order. More quantitatively, by expanding 
p2 — 1 — — a-iid performing the resulting Gaussian 
integrals, we find that 



s = 1 - 



1 



C2P2 ^1 



(37) 



to first order in \/ {p^l,n). To leading order, one has s = 1 
and Eq. (35) simplifies to 



h = Cp^ 



with 



C= 1 



C2 



Using this in Eq. (36) we get: 

f3 = lP^2 



where 



(6/5)ciC2Po(^^ 



C2 - Ci 



l+Cipo(P 



(38) 



(39) 



(40) 



(41) 



For our set of variables ^2^ '^^ ^'^^ have a set 

of three equations, consisting of the simplified Eqs. (38) 
and (40) and the osmotic pressure equality (26). Before 
turning to the latter, we note that Eq. (41) implies that 
P is of the same order as for — ^ 00 unless 7—^0 
in the same limit. We will find below that the second 
alternative holds, and that /3 is actually proportional to 

Consider now the pressure equality, given by Eq. (26) 
in its still exact form. We would like to apply to this the 
same procedure of splitting the /-integrals into contribu- 
tions from short and long rods, and then expanding in 
Pl and P2I in the short rods part. It turns out, however, 
that the expansion in the short rods part needs to be 
carried to third order in order to give an equation which 
is independent of the ones we have already derived. We 
therefore first rearrange the pressure equality into a form 
where the linear and quadratic terms in the expansion 
cancel automatically. 

On the l.h.s. of Eq. (26), we write po = Jdl p^°Hl); in 
the quadratic term we similarly replace one of the factors 
using Po — Pi ^ Jdl p^^\l)l. The r.h.s. can be rewrit- 
ten in a similar way, using the expressions (22,23,24) for 
the moments of the nematic phase. Including an overall 
factor of 2, the osmotic pressure equality then takes the 
form 

O^Jdlde p(o)(0 {e/5'+/(«)' [2 + cp^^l - C2P^IV2] 

- 2 - c^pol} (42) 

and the definitions (25) and (27) can be used to write the 
last two terms in the square bracket as 



cip^l - C2p^ir2 = [-f{0) -l3 + cipo]l 



In the same way, the definition (19) of /? can be written 
in integral form as 

P = jdlde p^"Hl)cil [1 - e^'+/(^)'] + C2P^ 

Multiplying by po — Jdl p^^^ {1)1 and bringing all terms 
to the r.h.s. gives 



0= fdlde p^°\l) [-pi + cpol [1 - e^'+^(«)' 



+ C2P^l} 

(44) 

If we now sum Eq. (42) and Eq. (44), and make use of 
Eq. (43) and the fact that C2P2 = — JdO f{9), we obtain 
after a few rearrangements 



0^ I dldd p^°\l)[ 

[P + f{0)] I [e^'+fW + 1] I 



(45) 



Eq. (45) is the desired form of the osmotic pressure equal- 
ity, which is still exact. We can now proceed as above 
and split the l-integral into the contributions from short 
and long rods. In the short rods part, we Taylor expand 
the integrand in the small variable [P + f{d)]l; as antic- 
ipated above, the first nonzero term in this expansion is 
of third order. To leading order, Eq. (45) then becomes 







dl dO p(^){l) [p + f{6)] 



6 



(46) 



where represents the long rods part of the integral in 
Eq. (45). In the first term we can again take the upper 
integration limit to 00 and obtain 



-^Po jde [p + f{e)r = -h 



(47) 



where as above the average (/'^) is over the normalized 
parent distribution P^^'>{1). 

Our final task in simplifying the pressure equality is 
to find I3 . Because pim diverges (and we have restricted 
the integration range to the long rods) the exponential 
factors in Eq. (45) dominate the integral. Of the terms 
involving these factors, the first one is also negligible com- 
pared to the second one, — [/? + f{9)]le^'''^^^^^'', since it 
is missing the large factor [P + f{0)]l. Comparing with 
Eqs. (32,33), we thus have that 

I3 = -[Ph+C2P^{l2 - h)] = -[P + C2p^is - l)]/l 

for large Eq. (37) further implies that C2p^(s — 1) = 
— l//„i to leading order, which is negligible compared to 
P since pim — ^ 00. Altogether we have the simple result 
I3 = —pii, and the transformed pressure equality (45) 
becomes in its final form 

-■^po fd9 [p + C2P^ {V2 - 1)] ' = Ph (48) 



(43) 
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C. The isotropic cloud point 

In the previous subsection, we used the large cutofF- 
limit to simplify the three basic equations for the three 
variables po, and /3, obtaining Eq. (38) for p2, Eq. (40) 
for (i and the transformed pressure equahty (48) for p^. 
While in principle these equations need to be solved si- 
multaneously, the limiting value of the isotropic cloud 
point is easy to extract directly, as follows. Using 
Eqs. (38,40) to eliminate [i and /i, and dividing by (p^) 
the pressure equality (48) becomes 

-^PoP? Jde [^ + C2 {V2 - l)f = C-1 (49) 

The l.h.s. of this equation tends to zero for l„i — *■ oo 
since p^ does; in the same limit we must therefore have 
C7 = 0. Thus either C = or 7 = 0; but the first case we 
have already excluded in our discussion after Eq. (31) as 
leading to the spinodal instability rather than the cloud 
point. Thus 7 = 0, and inserting ci = 2 and C2 = 5/4 
into Eq. (41) gives 



(50) 



the 



as our final result for the parent density po 
isotropic cloud point, in the large cutoff limit 00. 
From the definition (21) of the polydispersity a, and 
bearing in mind that we took the average rod length in 
the parent to be unity, we have {P) = cr^ + 1 so that the 
cloud point condition can also be written as 



1 1 



Po 



4cr^ + 1 



The result 7 — > for Z„, 00 implies, from Eq. (40), 
that P P2 for large Im- This relation helps to un- 
derstand the appearance of the minus sign on the l.h.s. 
of the transformed pressure equality (48). Call 17(6') = 
P + C2P2 0^2 — 1) the function appearing in the integrand. 
Then (3 = g{Q) > while for other 9 the value of g{9) is 
lower. In fact, if /? ^ p^, then g{9) is negative except in 
a very small region around 9 = (and of course 9 = tt). 
These negative values of g{9) will dominate the integral 
on the l.h.s. of Eq. (48) and cause it to be negative, mak- 
ing the l.h.s. positive overall. Interestingly, the situation 
on the r.h.s. is the reverse. As discussed above, for large 
Im we can write 

--/3 = fdldB p(°n0[/3 + /(6')]Ze^'+/(^)' 



dl d9 p^^\l)g{9)les^'^'^^ 



Since the integral here is over the long rods only, and in 
addition g{0)lm ~ plm ^ 00 for Im —> 00, the exponen- 
tial factor ensures that the integral is dominated entirely 
by the small region of angles around 9 = (and 9 = tt) 
where g{9) is positive. 



Having seen that <C p^ for large Im , we can now show 
that in fact (3 ~ (p^) . On the l.h.s. of Eq. (48) we can 
neglect the first term in square brackets and integrate the 

remainder exphcitly to obtain (54/35)((^'^)/6)poC2 (p2')'^- 
On the r.h.s. nothing changes and we can use Eq. (38) to 
eliminate Ii. Dividing by p^, Eq. (48) thus becomes 



showing that indeed /3 oc (p2')^ 



9_(p)po4 
35 C 



(51) 



D. Scaling of /3 and pa 

So far we have only used our simplified equations for 
the isotropic cloud point to deduce the value of the cloud 
point density for 00. The equations were obtained 

by expanding the short rod integrals in the equations 
for P2 and /3 to linear order in these two variables; in 
the osmotic pressure equality we needed to expand to 
third order because lower terms cancel out automatically. 
From the expanded pressure equality we deduced that 
the coefficient 7 relating /3 and p^ to linear order has to 
vanish in the limit — > 00, and that in fact /3 must be 
proportional to (p^)^ to leading order; the condition 7 = 
gave us the limiting value of the cloud point density for 
infinite cutoff. These results all apply generally for parent 
rod length distributions with fat (less than exponentially 
decaying) tails. 

We now proceed to obtain from our equations ad- 
ditional information regarding the variation of the key 
quantities with the cutoff Im, in the regime of large Im- 
These quantities are P, p^ and the moment densities of 
the isotropic cloud and nematic shadow phases. To do so, 
we have to explicitly calculate the long rod contribution 

h=po dl P''^\l)le('' /d^e'-^W 



to p5^, which so far we had left unevaluated. Anticipating 
that the integral is dominated by rod lengths of order Im, 
for which p^l is large since p^Un — ^ 00 as ^ 00, we 
can evaluate the angular integral by expanding around 
9 = (see before Eq. (37)) to get 



3C2P2 



(52) 



It will be useful to remove the overall factors and define 
the quantity 



J 3C2P^ 

ll — ll 

Po 



dl p(°)(/)e''' 



(53) 



which only depends on (3 and Z,„. 

v2 



From Eq. (38) we 

then have /i = (3c2C/po) (P2 and using Eq. (51) to 
eliminate p^ in favour of (3 gives 
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= A, 



3c2C 



(54) 



All factors in the definition of A, and therefore A it- 
self, tend to nonzero constants as /,„ — > oo; the above 
equation therefore implicitly determines the asymptotic 
dependence of /3 on Zm- 

To evaluate Ii, one notes (see App. B) that for fat- 
tailed parent distributions such as the log-normal the ex- 
ponentially increasing factor cxp(/3^) dominates the be- 
haviour of the integrand. The integral is therefore domi- 
nated by the region I k, and one can replace I = /,„ in 
the non-exponential factors of the integrand, yielding 



/i = p(°)(U- 



for large Ir, 



Eq. (54) then gives 



(55) 



(56) 



The function e^/x^ on the Ihs of Eq. (56) can be inverted 
for large argument x and yields asymptotically just a 
logarithm (see App. B), so that 



(57) 



Now for a log-normal parent distribution lnP*^'')(/) = 
— ln^^m/(2w^) plus terms that are negligible by com- 
parison for large Im, giving for the leading asymptotic 
dependence of (3 on 1^ 



In' L, 
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(58) 



For a more general parent distribution with a fat tail, i.e. 
P'^^^l) = e"''^'^ with h{l) a function that diverges with 
I less than linearly, the approximation (55) is not quite 
correct but can be refined (see App. B) to 



h 

Inserting into Eq. (54), 
obtains asymptotically 



fi-h'{l^) 



(59) 



one can again solve for (3 and 



hjlm) 



(60) 



which generalizes Eq. (58). For h{l) ^ I" with < a < 1, 
which corresponds to a parent distribution decaying more 
quickly than a log-normal but still less than exponen- 
tially, one gets P ~ ^m"^- Using that /3 oc (p^) for large 
Im, the above results also imply that ~ [h{lm)/lmV^^ , 
with P2 cx Irri^^'^ In/m for the log-normal case. Notice that 
the results for h{l) ~ P extrapolate sensibly towards the 
case of a parent distribution without a fat tail, i.e. a = 1. 



Then neither f3 nor have any reason to tend to zero for 
/,„ —^ 00, and this is consistent with the limit a — > 1 of 
the scalings above. As explained in Sec. Ill C, the scaling 
of P2 also determines directly the length-dependence of 
orientational order in the nematic phase: only rods with 
I around the "ordering length" and greater show 

significant ordering. For the log-normal case, the order- 
ing length is I/P2 ~ ^m^/ln^m and 3> 1 as expected, 
but simultaneously <C Im- This makes precise our ear- 
lier statements that orientational order is confined to the 
longer rods in the nematic phase; ordering appears for 
rod lengths far above unity, but still well below the cut- 
off length Im- 



E. Scaling of densities 

Finally, we study the scaling with of the isotropic 
cloud point density po and the number density p^ and 
rescalcd volume fraction pf of the nematic shadow phase. 
From Eqs. (29,30), with /? < and Iq <C h and insert- 
ing Eq. (38), one gets 



Po 
Po 



N 

C2P0P2 



[C ~ cpoil')] p^ 



(61) 
(62) 



where po is the density of the isotropic cloud phase at 
finite, large cutoff, which is still a function of Im- The 
interpretation of these results is that for finite but large 
Im, the cloud and shadow curves differ by terms of or- 
der P2 (~ Zm^^^ln/m for a log-normal parent), whether 
expressed in terms of po or pi. In the limit Im — ^ 00, 
the shadow curve collapses onto the cloud curve, and 
both tend asymptotically to the limiting form (50). This 
may appear counterintuitive, suggesting that the nematic 
phase does not at all "feel" the presence of the orien- 
tational order. The answer is of course that only the 
longer rods in the nematic are orientationally ordered, 
and for rod lengths comparable to the cutoff Im the ne- 
matic and isotropic density distributions do indeed differ. 
That there is a genuine difference between the isotropic 
cloud and nematic shadow phases becomes evident if we 
look at higher order moment densities / dl p{l)l" of the 
density distributions. In the parent, and therefore the 
isotropic cloud phase, all these moments will generically 
be finite (one exception being parent distributions with 
power-law tails). For the nematic phase, the picture is 
completely different. For n = 1 we have the moment p^ , 
and we showed above that the long rod contribution to 
this converges to zero for — > 00. For larger n > 1, since 
the long rod integrals are dominated by values of Z ~ ^ 
the long rod contribution will scale as l^^Ii ^ I] 
from Eq. (38) (see App. B). For the log-normal parent, 

we found that p^ ~ Im^^'^ Inlm and so all moments of 
the nematic shadow phase of order 71 > 3/2 diverge with 
Im- For parents with less fat tails, P*^°^ ^ exp(— Z"), we 

had instead f3 ~ l^-^ and p^ ~ (3^^^ - ^^n"^'^^ im- 
plying that nematic moments of order n > (3 — q;)/2 



m ; 
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diverge. Similarly, moment densities taking account of 
orientational order, / dl d9 p^{l, 9)l'^V2 diverge when the 
exponent n is large enough. As we can see then, the 
similarities between the isotropic cloud phase and the 
ncmatic shadow phase are rather more superficial than 
suggested by the fact that their low-order "physical" mo- 
ment densities po and pi agree, and that — > for 
— *■ oo; with respect to any quantity that attaches suf- 
ficient weight to the behaviour of the long rods, they are 
extremely different. 

To complete our analysis of the scaling of the moment 
densities, we need to know in Eqs. (61,62) how the cloud 
point density po depends on l^- By Taylor expanding 
the relation between 7 and po, Eq. (41), one easily sees 
that the deviation of po from its limiting value l/(4(^^)) 
is, to lowest order, proportional to 7. To determine how 
7 converges to zero, compare /3 = k {P2) from Eq. (51) 
with P = JP2 from Eq. (40). This shows that 7 should 
be proportional to p^, although we cannot determine the 
coefficient explicitly: the result /3 — 7p^ is only the low- 
est order in an expansion in p^ and we would need all 
terms contributing at quadratic order to get the propor- 
tionality coefficient in 7 oc p^. Nevertheless, 7 c>c p^ 
alone tells us that po approaches its limiting value as 



1 1 



Po 



4cr2 + 1 



-0(7) 



1 1 



4ct2 + 1 



+ 0{p^) (63) 



This result, together with Eqs. (61,62) shows that for 
large Im all the densities po, p^ and pf differ from their 
limiting value for Im — *■ cxd by terms proportional to p^. 



F. Numerical support for the theory 

The analytical results obtained so far are based on the 
fact that f3 and p^ must converge to zero as the cutoff 1^ 
becomes large; we supplemented this with the assump- 
tion that pim diverges, from which it then followed that 
also p^hn diverges. These statements allowed us to ar- 
gue that integrals over the rod length / could be split 
into two distinct parts: one for the short rods, where we 
could expand in the small parameters f3 and p^, and one 
for the long rods where the exponential factor exp{(3l) 
is large and dominant and leads to an increase in the 
nematic density distribution towards the cutoff length. 
We already saw some graphical support for this in Fig. 5 
but provide in this section more detailed numerical sup- 
port for the theory, focussing on the case of a log-normal 
parent distribution. 

From the theory we obtained that in the limit of in- 
finite cutoff, the cloud and shadow curves collapse onto 
the curve given by Eq. (50). In other words, the mo- 
ment densities po — pi of the isotropic cloud phase, and 
those of the nematic phase (p^ and p^) should tend 
to the same limit for /,„ — > 00. A plot of these den- 
sities against Im (Fig- 10), for fixed polydispersity cr, 
shows results compatible with this prediction, with all 



curves tending towards the theoretically predicted lim- 
iting value. A more detailed check is obtained by plot- 
ting the densities against p^. Our theory predicts that 
they should approach their common asymptotic value lin- 
early in p^, and Fig. 11 (a) confirms this. We can even 
check the known coefficients of this linear dependence 
by plotting the differences po — p^ and p^ — po against 
P2 ■ Eqs. (61,62) predict that po — Po — C2P0P2 and 
Pi — Po = [C — C2Po{i^)]P2 ■ In this representation, the 
l„i — > 00 limits of the coefficients of p^ on the right-hand 
sides are known, and the unknown proportionality factor 
in the approach of po to its limiting value (see Eq. (63)) 
has been eliminated. Figure 11 (b) displays the theoret- 
ically predicted linear relations and shows good conver- 
gence of the numerical results to these asymptotes. From 
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FIG. 10. Plot of the cloud point density po and the nematic 
shadow density pq and rescaled volume fraction against 
the cutoff Im, for fixed polydispersity a ~ 0.532 {w = 0.5). 
The horizontal line shows the common limiting value for 
Im — > 00 which our theory predicts, l/[4(cr^ + 1)] ~ 0.195. 



the theory we also obtained the scaling of /3 with giv- 
ing (3 ~ (In^ Im) / i'2lmw'^) to leading order. In Fig. 12 we 
check this by comparing plots of f3 and (In^ lj-n)/{2lmw'^), 
both against /,„. For a more detailed comparison, wc can 
estimate from our theory the first correction to the lead- 
ing behaviour of /?. From Eq. (57) wc see for a log-normal 
parent distribution that this correction term, which arises 
both from the factor ll^ in the denominator and the sub- 
leading term in lnP^^'{lm), should scale as ^ /^^ln/j„. 
(Again, because we have not systematically kept track of 
subleading terms, we do not attempt to determine the 
coefficient of this correction term.) This correction is 
smaller by 0{1/ \nlm) than the leading term, so that we 
should have 



PL 
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and a plot of f3lm/ ^'n^ Im against l/lnlm should show a 
linear approach to the limiting value l/(2u'^). The re- 
sults shown in the inset of Fig. 12 are compatible with 
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FIG. 11. (a) Same as Fig. 10, but with pa (which itself 
is a function of the cutoff Im) shown on the x-axis instead 
of Im itself. The results are compatible with the theoretically 
predicted linear approach, for p2 — > 0, to the common limiting 
value of 1/4((T^ + 1) « 0.195 (a « 0.532) that is marked by 
the dot. (b) Differences po — Po (bold solid) and Pi — po (bold 
dashed) versus p2 ■ For large Im (small p^), these converge 
to the linear relations predicted by Eqs. (61) and (62) and 
shown by the thin lines. 



this. The plot stiU shows some curvature, however, indi- 
cating that even our largest values of Im = 5000 are not 
yet large enough for higher-order corrections to be negli- 
gible. Accordingly, a linear regression gives an intercept 
of ~ 1.81, rather below the theoretical value l/2w'^ = 2 
for the chosen w = 0.5, while a cubic regression gives a 
much closer result of ~ 1.97. It is worth stressing that 
the above results support the assumption f3lm — *■ oo for 
Im oo which we made in our theory: the good fit to 
the theoretical asymptotics implies that f3lm ~ In^ Im for 
large Im, which indeed diverges. Even for smaller Im the 
trend for pim to grow with Im is noticeable; e.g. for the 
scenario shown in Fig. 12, (3lm increases from w 42 at 
= 50 to w 431 at Im = 5000. 

Finally, our theory also predicts how the parameter p^, 
giving the overall degree of orientational ordering in the 
nematic shadow phase, scales with Im- The result is that 
P2 should be proportional to /3^/^ for large Im, so that 
a plot of these two quantities against each other as Im is 
varied should give a straight line through the origin, with 
slope from Eq. (51). Fig. 13 is consistent with this, 
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FIG. 12. Plot of /3 (solid) and the theoretically predicted 
asymptotic scaling /3 — {In^ lm)/{2'w'^) (dashed) against Im at 
w — 0.5. The inset provides a more stringent test of the the- 
ory by plotting /3Zm/ln^?m versus l/lnZ™,. Theory predicts 
that this should approach l/2w^ — 2 (dotted line) linearly; 
the results plotted are broadly compatible with this but still 
show some curvature, indicating that our values of Im are 
still too small to have reached the asymptotic regime where 
higher-order corrections can be neglected. 



but shows the expected deviations in the pre-asymptotic 
regime where and (3 are not yet small. 



V. CONCLUSION 

We have analysed the effects of length polydispcrsity 
on the phase behaviour of hard rods, within the V2 On- 
sager model which is a simplified version of the full On- 
sager theory. We focussed in particular on rod length 
distributions with fat tails that decay less than expo- 
nentially, using log-normal distributions as our primary 
example. Our work was originally motivated by the ques- 
tion of whether, even for such unimodal length distribu- 
tions, the presence of relatively many long rods could 
cause I-N-N phase separation to occur. We do indeed 
find an I-N-N coexistence region in the phase diagram 
for a log-normal length distribution with a finite cutoff. 
Unusually, this region is located inside the I-N coexis- 
tence region, i.e. not connected to an N-N region, and 
very narrow in density, separating two different regimes 
of isotropic-nematic coexistence. At densities below the 
three-phase region, the nematic phase is significantly en- 
riched in long rods compared to the isotropic: its density 
distribution exhibits a second peak for large rod lengths, 
in addition to the one at lengths of order unity (the aver- 
age rod length in the system). Because the overall num- 
ber of long rods in the system is small, however, the ne- 
matic phase only occupies a very small part of the system 
volume in this regime. As the density increases towards 
the three-phase region, the second peak in the nematic 
length distribution moves towards the main peak; beyond 
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FIG. 13. Plot of pg* against yT?, obtained by varying Im at 
fixed w = 0.5 (solid), together with the asymptotic scaling 
P2 = {P/t^Y'^ predicted from Eq. (51) (dashed). As pre- 
dicted by theory, both (3 and tend towards zero as Im 
increases. In the same limit they also approach the predicted 
linear asymptote, while for smaller Im and therefore larger 
and P deviations from linearity occur. 



the thrcc-phasc region, the coexisting nematic only ex- 
hibits a single peak, and the phase behaviour resembles 
that found for systems without fat-tailed length distri- 
butions. Moving towards larger polydispersitics in the 
phase diagram, the three-phase region eventually closes 
off and is replaced by a crossover, within the I-N coexis- 
tence region, between the two types of nematic phase just 
described. Coming from small polydispersity ct, on the 
other hand, the three-phase region begins at a thresh- 
old value of (7 where the isotropic cloud curve shows a 
kink and the properties of the coexisting nematic shadow 
phase change discontinuously. 

In the theoretical part of the paper, we gave an ana- 
lytical treatment of the onset of nematic order, i.e. of the 
isotropic cloud point, in the limit of large length cutoffs 
Im. The key quantities for the theory were the param- 
eter (3, which determines the exponential enrichment of 
the nematic phase in the longer rods, and the average 
ordering of the nematic phase. We were able to show 
that both of these quantities must converge to zero for 
Im oo, and based on our numerical results supple- 
mented this with the assumption that pim — >■ oo. On 
this basis we could construct a self-consistent theory for 
the limit of large Im, by treating separately the behaviour 
of short rods (with lengths of order one) and long rods of 
order the cutoff length. Surprisingly, the theory predicts 
that the density (and rescaled rod volume fraction) po 
of the isotropic phase, and the density Pq and rescaled 
rod volume fraction of the nematic shadow phase, all 
converge to the same limiting curve when plotted against 
polydispersity a. This curve can be calculated explicitly 
and depends on no other details of the rod length distri- 
bution except for the presence of the fat tail. The theory 
in fact makes a stronger statement about the similarity 
of the isotropic cloud and nematic shadow phases: for 



any rod length of order unity, the density distributions 
in the two phases become equal for Im ^ oo, and the 
degree of orientational order of such rods in the nematic 
shadow tends to zero. This does not, however, imply that 
the isotropic and nematic phases coincide completely, or 
that the nematic lacks orientational order: when we con- 
sider the longest rods in the system, of the order of the 
cutoff length, we find that the nematic shadow contains 
a much larger number of them, which also exhibit strong 
nematic order. Our theory also allowed us ascertain the 
rate of convergence of all quantities to their limit values, 
to leading order as — *■ oo, and the predictions were in 
good agreement with numerical results. 

An important question that our current study cannot 
answer directly is how the above results for the V2 On- 
sagcr model will carry over to the full Onsager theory. 
The comparison of our earlier results for bidisperse dis- 
tributions'^*' were encouraging in this respect, showing 
qualitative agreement except for the presence of an N- 
N critical point in the 1^2 Onsager model which in the 
Onsager theory is replaced by an unbounded N-N coexis- 
tence region opening towards large densities. A direct at- 
tack on the full Onsager problem, with the required self- 
consistent determination of all orientation distributions, 
is of course a formidable numerical problem. Neverthe- 
less, we have recently obtained preliminary results which 
suggest that the isotropic cloud curve for Onsager theory 
with a log-normal length distribution indeed features a 
kink, as it does for the P2 Onsager model studied above; 
this shows indirectly that the full Onsager theory must 
also exhibit an I-N-N coexistence region as found above. 
The overall topology of the phase diagram remains an 
open question, however; in particular, one would like to 
know whether the I-N-N coexistence is again confined to 
a narrow interval or whether it extends over a broader 
range and is bordered by an N-N phase separation at 
higher densities. 
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APPENDIX A: THE CASE /? ^ 0, p^ FINITE 

In our theory for the isotropic cloud point we showed 
that in the limit of large cutoff Im ^ 00, the parameter 
(3 tends to zero while pim 00. The moment density 
P2 , on the other hand, which describes the orientational 
ordering of the nematic shadow phase, could in principle 
either also tend to zero, or converge to a nonzero con- 
stant in the limit. In this appendix we show that the 
second case can be excluded because it does not lead to 
a consistent solution of the cloud point equations. 

The assumption that f3 again allows us to break 
all integrals over I into two intervals, ^ = . . . ^ and I — 
I. . .Im with 1 <^ I Im- AH statements made in the 
main text regarding the long rod contributions lo , Ii and 
I2 to the nematic shadow phase moments Pq , p^ and p^ 
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also remain valid, seeing as the only property of that 
entered was that p^lm diverges as Im oo; if itself 
has a nonzero limit value then this is certainly the case. 
We know therefore that Iq is smaller by a factor 
than Ii and can be neglected, and that I2 = Ii to leading 
order since as a consequence of p^lm ^ 00 the longest 
rods are strongly ordered. 

The treatment of the short rod parts is actually easier 
than in the case of ^ 0: a finite limit value of 
ensures that these integrals have nonvanishing differences 
from the corresponding values in the isotropic phase, even 
for /„i — > 00. Any corrections from the small nonzero 
value of P are negligible by comparison, so that we can 
set /3 = directly and write Eqs. (22,23,24) as 



where 



N 7 



P0iJ2+l[) 



(Al) 



where I[ ~ h/pa and 



Ji^ Jdl dd /p(")(/)e^W' 
J2= fdlM lP^°\l)r2ef^'^^' 



Fip^) 



2(1- Jo) 
C2(Ji- J2 + I) 



It can be shown that this function passes through the 
origin with a slope of one, but then curves downwards 
and approaches a constant for large p^. It follows that 

Fip^) < P^ 

for p^ > 0, so that the only solution is p^ = 0, in con- 
tradiction to our assumption that p^ had a nonzero limit 
value for Im 00. This assumption must therefore be 
abandoned, and we conclude that p^ must tend to zero 
for Im ^ 00. From our arguments it is clear that this 
conclusion holds for any fat-tailed parent distribution. 

We mention finally that, even though the above treat- 
ment based on the assumption of nonzero p^ is not fully 
self-consistent, it actually recovers the correct limiting 
value of the cloud point density. As pointed out above, 
Eq. (A4) has p^ = as its only solution. Now for small 
P2 one can expand Jq = 1 — C2P^, Ji = 1 — C2P2 i^"^) 
J2 = (c2/5)p^(Z^); inserting these results into Eq. (A3) 
one finds po = 1/(4(Z^)) in the limit p^ 0, in agree- 
ment with Eq. (50). 



APPENDIX B: THE LONG ROD INTEGRALS 



(In principle, the short rod integrals run over / = . . . / , 
but because the integrals as written are convergent, one 
can take the upper limit to 00 without error at the level 
of our approximation.) We can determine the value of I[ 
from Eq. (27) for /3, which using P becomes 



= -ci (pf - Po) 
and solving for I[ one finds 

Cl — C1J1+ C2J2 



N 

C2P2 



Cl - C2 



(A2) 



Inserting this into the pressure equality (26), which with 
a common factor of po removed reads 

1 + yPo = Jo + yPo(Ji + I'lf - yPo(J2 + I[f 

we can solve explicitly to obtain the cloud point density 
Po as 



PO 



Cl - C2 2(1 - Jo) 

ClC2 l-(Jl-J2)2 



(A3) 



where the r.h.s. is now a function only of p^, via the 
short rod integrals Jo, Ji and J2. Finally, substituting 
the results (A2) for /( and (A3) for po into Eq. (Al) for 
P2 , we obtain after a little algebra the self-consistency 
equation 



P^ 



F{p^) 



(A4) 



In this appendix we discuss how to evaluate the general 
long rod integral 



dl P{l)re^^ / de e'/(^) 



Po 



with P{1) = P^"^l) for brevity. For n = and 1 this 
reduces to Iq and Ii, the long rod contributions to the 
moments Po and p^ of the ncmatic shadow phase. As 
explained before Eq. (52), since the integral will turn out 
to be dominated by the long rods with / « we can use 
the fact that these rods are strongly ordered {p^l ^ 1) 
to perform the angular integral and get 



Po 



3C2P2 



dl P{l)V- 



(Bl) 



Our task now is to evaluate /. If this really is domi- 
nated hy I K Im due to the presence of the exponential 
factor exp(/3^), we should be able to replace the rest of 
the integrand, f{l) = P(Z)/"~^ by its value at / = Im- To 
check the quality of this approximation, we include the 
first term in a Taylor expansion. 



/(O = f{lm) + lrn)r{U) + ■ 

which when inserted into Eq. (Bl) gives 



filn 



(B2) 



(B3) 
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If the leading term is a good approximation, the correc- 
tion term should be small, i.e. —f'{lm)/[l3f{lm)] — > 
for ^ cx). Whether this is true depends on the In- 
dependence of p. From Eq. (54), this is given implicitly 
by ///? = A, with A a constant and / evaluated at n — 1, 
i.e. for f{l) = P{1). Keeping the leading term of / only, 
we get 



(B4) 



Using the inverse of the function G{x) = jx^ on the 
l.h.s. this can be formally solved to give 



m 



Now from the definition of G{x) we have 

G{x) = y X = \ny + 2\nx 

and therefore by recursion 

X = \ny + 2 ln(ln y + 2 ln(ln y + 2 In . . .)) 

Asymptotically the first term is dominant, so that x = 
G~^{y) = \ny and the asymptotic behaviour of /? be- 
comes 



(3- 



(B5) 



The condition —f'{lm)/[Pf{lm)] for our leading 
order-approximation for / to be correct is thus 



film) 



/G™){lnA-ln[P(/„04]} 







(B6) 



If we now restrict our attention to parent distributions 
with finite mean, i.e. Jdl P{l)l < oo, then P{l„i)l'^^ 
for /,„ — > oo and the constant term In A becomes negligi- 
ble. We therefore need 



/,41n/)^(^„0 ^ /„[(lnPy(/,„) + (n-l)//, 
HP{lrn)ll] 



ln[P(/„0/2 



(B7) 



For a log-normal parent distribution, lnP(/) 
-(ln^/)/(2w2) for large / and so (InP)'(O 
— (In I) / {w'^ I) , and the condition becomes 



(In^ UV2 + 0(lnZ™) 



In L 







which is obviously fulfilled for — > oo. For power-law 
parents. In P{1) ~ \nl, one also finds that the ratio (B7) 
tends to zero logarithmically. In both these cases, one is 
therefore justified in only keeping the leading term in our 
Taylor expansion around I — Im, giving 



/ = 



13 



P{lrn)l 



This proves in particular the statement made in the main 
text that the effect of the weight function I" only comes 
in through the factor l^. It is also worth stressing that 
Eq. (B5) implies that f3l„i ~ — \n[P{lm)lfn] — *■ oo, con- 
sistent with our assumption of such a divergence. This 
result also justifies a posteriori the intuition that the in- 
tegrals for / and I arc dominated by the longest rods: 
the exponential integral giving the leading order approx- 
imation has significant contributions only where Im — I 
is of order so that the fractional deviation of the 
relevant Z-values from [Ij^ — l)/lm ~ ^/iPlm) tends to 
zero for Z,„ ^ oo. 

In the more general case of a less than exponentially 
decaying parent, say P{1) = e^''*^'^ and h{l) ~ I" with < 
Of < 1, we have (InP)'(O - and the condition (B7) 
becomes 



la I 



0{\nl^) 



a 7^ 



Here one therefore needs to go beyond the leading order 
approximation for /. To do so, one can write f(l) — 
exp[— /),(/) + {n~l) Inl] and keep in Eq. (Bl) linear terms 
in the Taylor expansion of the exponent around I = Im- 
This results in 

/ ^ g-hil^) + {n-l] lnl^+pi^ + {l~l^}[l3~h'{l^) + {n~l]/l^] 

in-l„-h{l,„)+IM^ 

For n = \ this gives Eq. (59) in the main text. (The 
term [n — l)/lm in the denominator can be neglected 
even for n ^ I since, having dealt with the power-law 
distributions case h{l) ^ Inl above, we can here assume 
that h{l) diverges more strongly than InL) From this one 
can again deduce the scaling of /3, using the condition 

i/p = A 



m-h'ilm)] 

To leading order this gives Eq. (60), (3 = h{lm)/lm] this 
can be shown by writing (3 = h{lm)/lm + substitut- 
ing into the above condition and verifying that the as- 
sumption that (5 is a small correction is self-consistent. 
Once the scaling of P is established, one can then show 
that higher-order corrections to the first-order approxi- 
mation (B8) are indeed vanishingly small as Im oo; we 
omit the details here. 
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In the following we will use the same convention as we 
used in the numerical solution of the phase coexistence 
equations, i.e., we will assume that the functional form of 
P^'^^l) is fixed as Im is varied. In this way the parent will 
be strictly normalized only in the limit for Im —> oo, and 
the average length will be strictly equal to 1 in the same 
limit. However, as we mentioned in Sec. Ill, the corrections 
to the strict normalization are small for sufficiently large 
cutoff. In the following, as we are interested in the large 
cutoff limit, we will therefore always assume pi = po for 
the parent distribution. 
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